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| Abstract. We report measurements of electric potentials at the surface of a spherical 

£NJ ■ container of liquid sodium in which a magnetized inner core is differentially rotating. The 

azimuthal angular velocities inferred from these potentials reveal a strong super-rotation of 
' the liquid sodium in the equatorial region, for small differential rotation. Super-rotation was 

observed in numerical simulations by Dormy et al. . We find that the latitudinal variation 
of the electric potentials in our experiments differs markedly from the predictions of a similar 
numerical model, suggesting that some of the assumptions used in the model - steadiness, 
equatorial symmetry, and linear treatment for the evolution of both the magnetic and 
velocity fields - are violated in the experiments. In addition, radial velocity measurements, 
using ultrasonic Doppler velocimetry, provide evidence of oscillatory motion near the outer 
sphere at low latitude: it is viewed as the signature of an instability of the super-rotating 
O ■ region. 

' Keywords. Spherical Couette flow, super-rotation, magnetostrophic, liquid sodium 

!> ■ experiment, Hartmann layers, ultrasonic Doppler velocimetry. 

I> ' 

1. Introduction 

' In the recent years a revived interest has been shown for the rotating spherical Couette flow: 

the flow between two concentric spheres that rotate around a common axis z with different 
angular velocities ft (outer sphere) and + AO (inner sphere). Among the reasons for this 
• i— i ' renewed interest are the possibility of super-rotation when a magnetic field is added, and 

the prospect of observing dynamo action with such a flow. 

The first studies [2] dealt with the axisymmetric flow in the limit Ro = AQ/fl — ► 
and E = vjVLa? — > 0, where Ro and E are the Rossby and Ekman numbers, respectively, 
v is the kinematic viscosity of the fluid, and a the radius of the outer sphere. Proudman 
^ . P] discovered that, under the action of the Coriolis force, a vertical shear layer sets in to 

accommodate the transition between the fluid at rest (in the frame rotating at O) outside 
the vertical cylinder tangent to the inner sphere, and the fluid rotating at an angular velocity 
ranging from + AO, to f2 + Afi/2 inside the tangent cylinder. Stewartson [2] unraveled the 
inner structure of this shear layer for small yet finite values of the Ekman number, showing 
that an ageostrophic layer of width aE 1 / 3 is nested between two geostrophic layers of width 
aE 2 / 7 (inwards) and aE 1//4 (outwards). Years later, the validity domain of these analytical 
results was established in numerical simulations [U ^ . 

The Stewartson shear layer is expected to become unstable to non-axisymmetric per- 
turbations when the Rossby number is increased above a critical value Ro c . Experimental 
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evidence was found by Hide and Titman j^j- These instabilities were also studied in 3D- 
numerical simulations by Hollerbach who emphasized the difference between the Ro > 
and Ro < cases. When the Ekman number is small enough then Ro c <C 1, and experi- 
ments demonstrate that the instabilities are ^-independent |Sj . Schaeffer and Cardin [S] 
thus conducted a quasi-geostrophic analysis, and showed that the instabilities are akin to 
Rossby waves, thereby explaining the dependence on the sign of Ro. They found that the 
onset scales as: Ro c oc /3E 1 / 2 , where (3 is associated to the slope of the sphere at the lati- 
tude of the Stewartson layer. Above the threshold, a strong anisotropic and inhomogeneous 
Rossby wave turbulence develops, which exhibits a rather steep m -5 spectrum of kinetic 
energy in the azimuthal direction [S] (m is the azimuthal number). It means that a direct 
cascade of energy is prohibited by the geostrophic constraint, leading to a very low kinetic 
energy of the fluid flow at small scales. 

Once the instabilities are strong enough and when the fluid is an electric conductor, 
the instabilities can produce a magnetic field through dynamo action This happens 

when the magnetic Reynolds number Rm = Uajr\ is larger than some critical value Rm c 
(where U is the typical flow velocity of the fluid and r\ its magnetic diffusivity) . Within 
the quasi-geostrophic approach, Rm c is found to be in the range of 3000. However, full 
3D numerical simulations yield much smaller values around 100 |llj . but they are limited 
to larger values of both the Ekman number and the magnetic Prandtl number Pm = v/rj. 
These results validate the project of building an experimental dynamo based on the rotating 
spherical Couette flow • After the success of the liquid sodium dynamos in Riga an d 
Karlsruhe ^3], such a flow would get us one step closer to situations relevant to planetary 
cores. 

However, the presence of a magnetic field B will have a strong influence on the flow 
itself. This is well illustrated by the numerical studies of Hollerbach jlj and Dormy et al. 
PQ, who investigate rotating spherical Couette flow with an imposed magnetic field. Both 
studies only deal with the axisymmetric solution obtained in the limit Ro — > ((u ■ V)u 
terms neglected in the Navier-Stokes equation) and Rm — > (effects of the induced magnetic 
field are neglected). One effect of the magnetic field is to thicken the Stewartson layer. For 
large enough magnetic fields, the quasi-geostrophic character of the flow is broken. More 
intriguing, Dormy et al. [y discovered that, when the imposed magnetic field is an axial 
dipole, the fluid can rotate at an angular velocity larger than that of the inner sphere, in 
particular when the outer sphere is at rest. This phenomenon of super-rotation is illustrated 
in Figure n which shows the dimensionless angular velocities (CT^/sAO, where U v is the 
azimuthal velocity and s the cylindrical radius), meridional electric currents and electric 
potentials in a meridional plane. A wide crescent of super-rotation is clearly visible in 
Figure This phenomenon may be explained using the circulation of electric currents 
(Figure |3)) : the differential rotation generates a viscous shear at the outer boundary and 
electric currents are induced in the shear layer (named Hartmann layer) so that Lorentz 
forces balance the viscous forces. The electric currents of the Hartmann layers loop into 
the bulk of the fluid. They tend to follow the dipolar magnetic field lines but this is not 
possible in the equatorial region where they have to cross the magnetic field lines. In this 
region, they generate Lorentz forces that accelerate to a super rotating angular velocity a 
ring of fluid around the magnetic field line tangent to the outer sphere. An asymptotic 
analysis of this peculiar phenomenon has been carried out by Dormy et al. ^S] for $1 = 
(also see an d EZj)- R indicates that the width of the detached shear layer varies as 
Ha -1 / 2 , while the thin layers on the rigid surfaces are classical Hartmann layers [2B] whose 
thickness varies as Ha" 1 (the square of the Hartmann number Ha characterizes the ratio of 
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Figure 1: Super-rotation in axisymmetric numerical simulations (outer sphere at rest). The 
numerical method is described in ^H] and has been extended to take into account the finite 
electrical conductivities of the outer shell and inner sphere of the DTS experiment. The 
strength of the magnetic field (Hartmann number) is also the same as in DTS (see Table 
|5J). Isovalue contours in a meridional plane of (a) the normalized angular velocity U^/sAQ, 
(maximun 1.53, contour level 0.1, dashed line for value lower than 1, dotted line for value 
1 and continous line for value larger than 1); (b) the dimensionless stream function of 
meridional electric currents (maximun 0.93, contour level 0.1) - the electric currents are 
clockwise; (c) the normalized electric potential V/a 2 AflB in the frame of the outer sphere 
(maximum 600, contour interval 30). 

electromagnetic forces to viscous forces, and is defined in table |2j). 

Hollerbach |18U19j showed that the resulting flow and the phenomenon of super-rotation 
strongly depend upon both the geometry of the applied magnetic field (axial versus dipolar 
for example) and the electric boundary conditions on the inner and outer spheres. In 
particular, when both boundaries are perfectly conductive and for an axial homogeneous 
magnetic field, Hollerbach JH] and Hollerbach and Skinner [201 observe that the fluid rotates 
in the direction opposite to that of the driving inner sphere, and that the amplitude of this 
counter-rotation increases with the strength of the applied magnetic field. Buehler 
provided a simple explanation for this phenomenon: he showed that the flow rate in the 
counter-rotating jet is independent of the Hartmann number (provided it is large enough), 
and since the width of the jet decreases as Ha -1 / 2 , its velocity must increase as Ha 1 / 2 . 

In contrast, with a dipolar magnetic field and an insulating outer sphere, the super- 
rotation is positive and remains below 1.4 for £1 = [TJ. Hollerbach and Skinner [2U] also 
studied the linear stability of the axisymmetric jets for an applied axial magnetic field. 
Whatever the electrical boundary conditions, they predicted the onset of instability for 
Reynolds numbers of a few thousands and their non-linear calculations reveal a few drifting 
vortices aligned with the imposed magnetic field. 

When global rotation is present ($7^0), a competition takes place between the Lorentz 
and Coriolis forces. As the Elsasser number A = aB 2 / pQ decreases (a the electrical con- 
ductivity and p the density of the fluid), the ring of super-rotation stretches in the vertical 
direction and gets closer to the tangent cylinder, while super-rotation decreases The 
magnetostrophic regime, in which both forces are comparable is of particular interest, since 
it is the regime expected for planetary dynamos (e.g. |22j). The above analyses always 
consider very low magnetic Reynolds numbers where the non-linear action of the induced 
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magnetic field is neglected. In an attempt to retain non-linear contributions in the Lorentz 
force (Rm 7^ 0), while neglecting those in the advection term (Ro — ► 0), Cupal et al. [23] 
found a substantial increase of super-rotation up to 3.5 in numerical simulations. 

Little is known when all non-linear effects are present. Is super-rotation still present? 
How is the angular velocity distributed? Which symmetries are preserved? What are the 
characteristics of the meridional circulation? 

In this paper, we present the first experimental measurements of velocities in a magne- 
tostrophic spherical Couette flow in presence of a strong dipolar magnetic field. Azimuthal 
velocities in the liquid sodium flow are obtained by measuring the electric potential at the 
outer surface. Radial profiles of the radial velocity are measured with an ultrasonic Doppler 
velocimeter. A wide range of parameters is explored, with magnetic Reynolds number Rm 
up to 25. Experiments in a similar geometry have been conducted with a weak axial mag- 
netic field and without global rotation |24l I25j . Turbulent liquid sodium flows in cylinders 
have also been investigated [2T>] . 

The DTS set-up and measurement techniques are described in section 2. The relevant 
dimensionless numbers are listed in tabled Section 3 is devoted to our experimental results, 
which are discussed in section 4. 

2. Experimental set-up 

A schematic drawing of the central part of the DTS experiment is shown in Figure El 
Forty litres of liquid sodium are contained between a 7 cm-radius inner sphere and a 21 
cm-radius outer spherical surface. The 5 mm-thick outer spherical shell is made of stainless 
steel, which is about 8 times less electrically conductive than sodium at 130° C, the typical 
experimental temperature. The relevant properties of liquid sodium at 130°C are given in 
Tabled The inner sphere is made of two hollow copper hemispheres filled with 5 layers of 
rare-earth cobalt bricks. The copper wall thickness varies from 10 to 20 mm. The electrical 
conductivity of copper is about 4 times larger than that of sodium at 130° C. Once filled, the 
two hemispheres have been welded using an electron beam technique at the Laboratoire de 
l'Accelerateur Lineaire in Orsay. The assembly has finally been magnetized in a 5 T axial 
magnetic field, at the High Magnetic Field Laboratory in Grenoble. Outside the copper 
sphere, the magnetic field is found to be an axial dipole, aligned with the axis of rotation 
and pointing downwards, with: 

B(r,0) = ^-(2e r cos9 + e sm9) 
47rr 

in spherical coordinates, with 9 the colatitude and M = —700 Am 2 the magnetic dipolar 
moment, yielding a magnetic field magnitude ranging from 0.41 T at the poles of the inner 
sphere down to 0.008 T at the equator of the outer sphere. 

2.1 Electric potentials and azimuthal velocity 

Electric potentials are measured at several latitudes on the rotating outer sphere (see Figure 
|2|. The electric potential is obtained by simply sticking copper electrodes into 1 mm- 
diameter 4 mm-deep blind holes drilled in the stainless steel wall. Stainless steel is a poor 
electric conductor as compared to liquid sodium and thus it does not affect the dynamically 
generated electric potential, but it is a good enough electric conductor for the impedance 
of electrode pairs to be much smaller than that of the acquisition system. The signals pass 
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Figure 2: Sketch of the core of the DTS experiment. The 7 cm-radius inner sphere is made 
of copper and contains 5 layers of magnets. The inner radius of the stainless steel outer 
shell is 21 cm and its thickness is 5 mm. Tiny holes are drilled along several meridians and 
parallels to receive copper electrodes. They are 1 mm-diameter and 4 mm-deep blind holes 
as shown in the blown-up cross-section on the right. The large holes at latitudes —20°, 10° 
and ±40° receive interchangeable assemblies; a cross-section through one of them is shown 
on the left. An ultrasonic transducer can be mounted inside it. 



symbol property units value 



p 


density 


kg m 3 


930 


a 


electrical conductivity 


Q- 1 m- 1 


9 10 6 


V 


kinematic viscosity 


i —1 
m s 


6.5 1(T 7 


V 


magnetic diffusivity 


2 —1 

m s 


0.087 


c 


sound velocity 


m s^ 1 


2600 



Table 1: Physical properties of liquid sodium at 130° C (from |27|). The magnetic Prandtl 
number of sodium is Pm = u/r] = 7.5 10~ 6 . 
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through slip-rings for acquisition in the laboratory frame, where they go through a simple 
RC anti-aliasing 215 Hz low-pass filter, before being sampled at 1000 sample/second with 
a PXI-6229 National Instruments acquisition card. 

The electric potentials measured at the surface can be related to the azimuthal velocities 
beneath the Hartmann layer. Indeed, from the ^-component of Ohm's law: 

j = a(E + UxB) (1) 

we deduce (assuming steady-state): 

* < 2 > 

where AV is the electric potential difference between two electrodes placed along the same 
meridian and separated by a latitudinal angle AO, XJ^ is the azimuthal velocity of the fluid, 
B r the radial component of the magnetic field at the latitude of the electrodes, and je the 
^-component of the electric current. Further assuming that the latitudinal component of 
electric currents is negligible outside the Hartmann layer, and since the electric potential is 
constant across Hartmann layers |28j . the potential difference measured at the surface can 
be directly related to the bulk velocity: 



AV = a AO B r Up (3) 

Note that with our sign convention and since the magnetic dipole has its south pole in 
the northern hemisphere, we expect positive AV in that hemisphere when the inner sphere 
spins from West to East (U^ > 0). 



2.2 Ultrasonic Doppler velocimetry and meridional circulation 

We measure radial profiles of the radial velocity of the fluid by ultrasonic Doppler velocime- 
try. At two different latitudes (10° and —20°) on the rotating outer sphere, we install 
ultrasonic transducers in removable assemblies (see Figure EJ). The transducers operate at 
4 MHz and their active part is 5 mm in diameter (type TR0405LS from Signal Processing, 
Switzerland). The slightly diverging ultrasonic beam enters liquid sodium after passing 
through a 1.4 mm-thick stainless steel wall, which is polished on the sodium side. 

The electric signals pass through the above-mentioned slip-ring system, to be analyzed 
in the laboratory frame. Note that the 36-track system is enclosed in an electromagnetic 
shield, to prevent contamination of the weak back-scattered echoes by ambient noise. The 
treatment of the ultrasonic pulses and echoes is performed on a DOP 2000 velocimeter 
(Signal Processing). The instrument issues a series of 8-cycles pulses and listens to echoes, 
which are back scattered by weak heterogeneities present in the fluid. The time elapsed 
between the pulse and the echo provides the distance of the scatterer, while the associated 
Doppler frequency shift yields the radial velocity at that distance. In this manner, radial 
profiles of the radial component of velocity can be obtained, and repeated every 30 ms or 
so. Since liquid metals are opaque, ultrasound Doppler velocimetry is the most promising 
non-intrusive technique for mapping flows in MHD and dynamo experiments. However, the 
application to liquid metals raises specific issues and proves difficult: successful measure- 
ments in liquid gallium and sodium are recent [2^11301 ■ While electric potentials at the 
surface of the sphere are sensitive only to the azimuthal velocities, we get constraints on 
the meridional circulation from the ultrasonic Doppler velocimetry. 
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symbol 


expression 


units 


value 




a, 


outer radius 


in 


0.21 




b 


inner radius 


in 


0.07 




T l 

J 


a 2 /ir 2 r] 


s 


0.05 










B = Bi 


B = B n 


Ha 


aB 1 ^/fiopur) 


m s _1 


5200 


210 


U n 

*- ft 


B 1 ,/LlnO 


5.8 


0.2 




a/ ' UnPvri 1 B 


in 


4 10" 5 


10~ 3 








/ = 1Hz 


/ = 5Hz 


E 


u/Q.a 2 




2 10~ 6 


5 10" 7 




\/u/Q 


in 


3 10~ 4 


1.4 10~ 4 








B — B{ B — B B 


= Bi B = Bq 


A 


aB 2 /pQ, 




62 0.1 


12 0.02 


A 


V„/att 




4.4 0.2 


0.9 0.04 








A/ = 1Hz 


Af = 5Hz 


U 




m s _1 


0.44 


2.2 


Rm 


U a/r] 




1.1 


5.3 


Re 


Ua/u 




1.4 10 5 


7 10 5 








B = B t B = B B 


= Bi B = B Q 


N 


aaB 2 /pU 




180 0.3 


37 0.06 



Table 2: Typical values of the relevant parameters and dimensionless numbers for different 
imposed rotation frequencies / = Q/2-7T and differential rotation frequencies Af = AQ/2tt. 
For the numbers that depend on the magnetic field strength, two values are given, the 
first one with B = Bi = 0.2 T at the equator of the inner sphere, the second one with 
B = B a = 0.008 T at the equator of the outer sphere. 

2.3 Dimensionless numbers 

Table El gives the typical values of the dimensionless numbers in our experiments. We note 
that the Joule dissipation time tj is very short. The Hartmann number Ha does not depend 
on the imposed rotation rates: it is large everywhere in the sphere, and the thickness of 
the regular Hartmann layers <5na is less than 1mm. U a is the velocity of the Alfven waves. 
When the outer sphere spins, the Ekman number is less than 2 10 -6 and the thickness of 
the regular Ekman layers less than 0.3 mm. The values of the Elsasser number A and of A 
indicate that the magnetostrophic regime (A ~ 1 and A <S 1 ^2]) can be achieved in the 
interior of the sphere, while the flow is strongly influenced by the magnetic field near the 
inner sphere, and by global rotation near the outer sphere. Even for moderate differential 
rotation (Af = 1 Hz), the magnetic Reynolds number Rm is larger than 1, meaning that 
the effective magnetic field depends on the flow itself. Note that the Rossby number Ro is 
simply defined as Ro = = Af/f. Finally, the values of the interaction parameter 

N indicate that the flow is strongly influenced by the Lorentz forces, except near the outer 
sphere, where (u ■ V)u terms dominate. However, we note in the discussion section that 
even with a strong magnetic field, some inertial terms can play an important role. 



7 



> 
E 

CD 
O 

c 

CD 

o 
3= 



re 
+■« 
c 
o 

o 

Q. 

U 



5 -1 



> 

E, 

o 
u 
c 
o 

o 
3= 



.2 

c 

o 

o 
a 

a 




30 



20 



10 



N 



■10 



-20 







1000 2000 
time (s) 



3000 



-1 




1000 2000 
time (s) 



3000 



Figure 3: Temporal records of the differences of electric potential between pairs of electrodes 
placed at different latitudes at the surface of the outer sphere. The rotation rate A/ of the 
inner sphere is raised in steps (dashed lines and right hand-side vertical axis), (a) outer 
sphere at rest (/ = 0). (b) rotating outer sphere (/ = 4.5 Hz). 
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3. Results 



3.1 Azimuthal velocities from electric potentials at the surface of the 
outer sphere 

Figure |31 shows temporal records of the differences of electric potential between several pairs 
of electrodes at the surface of the outer sphere. In Figure I3K the outer sphere is at rest, 
and the rotation rate A/ of the inner sphere is raised in steps between -8 Hz and 8 Hz. 
During each step, the rotation rate is kept constant to within 5 10~ 3 during some 300 s. 
The electric signals are low-pass filtered to retain frequencies below 10 Hz. Each velocity 
step generates a jump in electric potential, forming a plateau sometimes preceded by a 
slow relaxation (probably of thermal origin). The signals in Figure |3] demonstrate that the 
fluctuations of azimuthal velocity in this frequency range are very small compared to their 
mean value, and that the coupling between the inner sphere and the fluid is very strong 
(between steps, the fluid reaches its equilibrium velocity in less than 2 s). The amplitude 
of the signals is in the mV range. 

In each pair, the electrodes are all 10° apart, and the latitude of their mid-point ranges 
from 10° to 40°. For 40°, two pairs have been installed along two different meridians. 
Their signals are undistinguishable in Figure 131 The records show that the difference in 
electric potential increases with latitude, related to variations of the velocity field as the 
corresponding increase is known for the imposed radial magnetic field B r (see equation 0J. 

Figure |3t> shows a similar temporal record, but this time the outer sphere spins at a 
constant rotation rate of 4.5 Hz during the whole sequence. The same behavior is observed. 
For each imposed A/ we retrieve the value of the electric potential plateau. 

We first concentrate on the signal measured at latitude 40°. We derive U v , using equa- 
tion |31 and divide it by sAO , where s is the cylindrical radius at that latitude. We can thus 
compare the azimuthal velocity of the fluid relative to that of solid body rotation (with the 
inner sphere) for different / and A/. 

Figure 0] summarizes our results as a function of A/, for various /. We observe that, for 
a given /, the normalized velocity is maximum, reaching 0.8 to 0.9, when A/ is minimum. 
We note that azimuthal velocity levels off at a higher value when global rotation is present, 
in particular for Ro > 0. 

We now discuss the case without global rotation (/ = 0). The normalized azimuthal 
velocity decreases rapidly as |A/| increases, and seems to level off at a value of 0.5 for 
|A/| > 10 Hz. Our observations can be extrapolated to the Ro — > (A/ — > 0) limit, for 
which numerical and analytical results are available. 

The numerical simulations provide an excellent guide to interpret our results. We have 
run the numerical model described in JH]) with parameters as close as possible to the 
experimental conditions: same Hartmann number (see Table |2J), electrical conductivities in 
the ratio 32/8/1 for inner sphere/fluid/outer shell, with appropriate thickness for the latter. 
The geometry of the applied magnetic field is the same as in our experiments (axial dipole). 
Figure^shows the contours of U^/sAQ, meridional electric currents, and electric potential 
V obtained numerically. At latitude 40°, the Hartmann layer is well identified and is very 
thin. The actual azimuthal velocity beneath that layer at 40° is given for comparison with 
our experimental results in Figure0]('*' symbol). 

We now investigate how the electric potential differences vary with latitude. The first 
move is to simply get U v /sA£l for the other latitudes using equation [31 This is what we 
have done in Figure for the experiments with the lowest |A/| and / = 0. The data 
show that the angular velocity of the fluid increases towards the equator, reaching super- 
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Figure 4: Variation of the azimuthal velocity (normalized to that of solid body rotation 
with the inner sphere) at latitude 40°, as a function of the rotation rate Af of the inner 
sphere, for various rotation rates / of the outer sphere. Two independent series (1) and 
(2) are presented. Typical uncertainties are below 10% within a given series, but different 
series can be separated by a slightly larger amount. The star (model) at Af = is derived 
from the numerical model for / = as discussed in the text. 
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20 40 60 80 

latitude (°) 

Figure 5: Comparison between the azimuthal velocity beneath the Hartmann layer (solid 
line, r/a = 0.95) and the azimuthal velocity inferred from the surface electric potential 
(dotted line, r/a = 1) using equation |21 as a function of latitude, for the same numerical 
simulation as in Figure ^ The dashed line is the velocity we would derive from the electric 
potential at r/a = 0.95. The symbols with error bars are the experimental normalized 
angular velocities, inferred from the surface electric potential, for the lowest |A/| with 

/ = o. 
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Figure 6: Latitudinal profiles of the differences of electric potential relative to that at 
latitude 40° for various rotation rates A/ of the inner sphere, (a) outer sphere at rest 
(/ = 0). The solid line is the prediction of our numerical model, (b) rotating outer sphere 
(/ = 4.5 Hz). IThe dashed line indicates the potentials that would result from global solid 
body rotation with the angular velocity observed at latitude 40°. 

rotation angular velocities up to 30% larger than solid body rotation with the inner sphere. 
However, the assumptions used to derive equation |31 appear to break down when approaching 
the equator. This is illustrated by the curves derived from the above-mentionned numerical 
model. The true azimuthal velocity at r/a = 0.95, outside the Hartmann layer, exhibits a 
slow decrease from the high latitudes down to 20°, where it rises sharply to a maximum of 
more than 1.3 at the equator. In contrast, the azimuthal velocity inferred from the electric 
potential at the surface using equation |3] is almost twice smaller at the equator, while the 
differences are small above 40°. There are two causes for this discrepancy and the clue is 
given by the third curve, which shows the azimuthal velocity that would be derived from 
equation |21 using the electric potential at r/a = 0.95. This time, the agreement with the 
true velocity is perfect down to latitude 20°. Around the equator, the velocity is again 
underestimated. This demonstrates that the contribution of electric currents in the core 
(see equation I2J) , is indeed negligible at latitudes above 20°, but that the electric potential 
is not quite constant across the Hartmann layer. In the Tropics, both approximations break 
down: the electric potential is not constant across the Hartmann layer, and the electric 
currents are not negligible beneath it. It is clear from Figure ^ that we are still under the 
influence of the Ha -1 / 2 thick layer at these latitudes. 

Going back to our experimental measurements, we thus decide to directly plot the elec- 
tric potential differences as a function of latitude. However, in order to focus on latitudinal 
variations, we normalize them with the potential difference measured at 40°. Figure 
shows our results for / = and A/ varying from —8 to 8 Hz. All profiles are remarkably 
similar despite a factor of almost 2 in their absolute level (see Figure |1} . The line derived 
from our numerical model falls well below our experimental points. 

The measurements for / = 4.5 Hz are given in Figure|Ht>. A similar behavior is observed. 
Only those points with |Ro| = \Af /f \ > 2 depart from the general trend, displaying smaller 
relative electric potential differences near the equator, below the line that corresponds to 
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solid body rotation with the angular velocity observed at latitude 40°, according to equation 

01 



3.2 Meridional circulation from ultrasonic Doppler velocimetry 

We now turn to the constraints brought on the meridional circulation by the analysis of 
radial profiles of radial velocity, obtained by ultrasonic Doppler velocimetry. Note that we 
cannot measure azimuthal velocities with this technique in the present set-up. Figure 
shows two examples of spatio-temporal plots and corresponding probability density func- 
tions (pdf) of velocity as a function of depth. Two different behaviors are observed. When 
the outer sphere is a rest (top row), strong and rapid oscillations extend from the outer 
surface down to about 40-50 mm. The corresponding pdf is very asymmetric: centripetal 
velocities reaching much larger values than centrifugal ones. Smooth stable centripetal ve- 
locities are measured from 80 mm down to the inner sphere. The equatorial oscillations 
are reminiscent of the jet instabilities observed in spherical Couette flow in the absence of 
a magnetic field. The rotation of the outer sphere kills the equatorial oscillations (bottom 
row), and the pdf is Gaussian at all depths. 

Getting back to the case with an outer sphere at rest, we find that the amplitude of the 
equatorial oscillations increases with |A/|. Around the depth of 30 mm, we measure the 
spread AU r for which the pdf has dropped down to one tenth of its maximum value. When 
the pdf are rescaled using AU r and an offset Uq, they all have the same shape, as shown in 
Figure |S^. We then plot AU r /sAS} as a function of Af (Figure [SJd) and observe the same 
kind of pattern as for azimuthal velocities: the ratio is maximum for small values of Af. 
The ultrasonic probe has been installed at two different latitudes: +10° and —20°. When we 
compare the two series of results, we find that higher AU r /sQ values characterize latitude 
+10°. Since the equatorial oscillations probably result from instability of the equatorial 
azimuthal jet, this observation is in agreement with an increase of differential rotation 
towards the equator. 

When the outer sphere rotates, the equatorial oscillations disappear, even when the 
Rossby number Ro = Af j f gets down to -2. 

4. Discussion 

We have conducted experiments on the spherical Couette flow with global rotation and an 
imposed dipolar magnetic field. The working fluid is liquid sodium, and we have developed 
several instrumental techniques to characterize the flow. Surface electric potentials convey 
crucial informations about the azimuthal velocities in the fluid. At latitude 40°, the di- 
mensionless angular velocities Um/sAQ we retrieve reach 0.9. They are maximum for small 
differential rotation rates. When the outer sphere is at rest, the extrapolation to Af = is 
in excellent agreement with the predictions of our numerical model. This numerical model 
exhibits super-rotation. It uses the same dipolar magnetic field and Hartmann number as 
in the experiments, and takes into account the finite electrical conductivities of the inner 
and outer walls. 

However, the latitudinal variation of electric potentials measured in the experiments 
differs markedly from the numerical predictions, and exhibit larger-than-expected values at 
low latitudes, which are not explained by the finite conductivity of the outer wall (correctly 
taken into account in the numerical model). Under the assumptions of the numerical model 
- steadiness, equatorial symmetry, absence of non linearity for the magnetic field (Rm <C 1) 
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Figure 7: Ultrasonic Doppler velocimetry results: spatio-temporal plots and histograms of 
radial velocity as a function of depth. The ultrasonic transducer is at latitude —20° and at 
depth=0 at the top of the plot. The ultrasonic beam hits the inner sphere at a depth of 
139 mm. On the left, the grey tones give the amplitude of radial velocity (scale in mm/s), 
centripetal in white and centrifugal in black. On the right, the grey tones are levels of 
the probability density function (pdf) of velocity, (a) and (b) outer sphere at rest and 
A/ = 1.2 Hz. Note the strong oscillations around depth=30 mm. (c) and (d) / = 4.7 Hz 
and A/ = 3.0 Hz. 
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Figure 8: (a) normalized pdf of radial velocities for / = at a depth of 27 mm, from 
ultrasonic Doppler velocimetry. Once rescaled with each width AU r and shift Uq, all pdf 
collapse on a single curve, (b) AU r /sAQ as a function of A/. Circles: / = 0. The full 
symbols are for a transducer at a latitude of —20°, while it is at 10° for the open symbols. 
The arrows indicate the shift to be applied to the data because different instrumental 
settings have been used for these data points. The squares are for / = 4.5 Hz and show 
that the equatorial oscillations are suppressed by global rotation. 

and for the velocity field (Re < 1) - these electric potential measurements indicate a much 
stronger super-rotation of the fluid than predicted numerically. 

The difference could be due to the meridional flow: while it is absent from the numerical 
solution for / = 0, we measure radial velocities reaching about 10% of the azimuthal ve- 
locities. The meridional flow advects the magnetic field lines, producing additional Lorentz 
forces. Thus, the meridional flow could be driven by centrifugal forces (one component of 
the Reynolds stresses) resisted by magnetic forces. According to this scenario, magnetic 
field lines and the associated region of super-rotation would be pushed closer to the surface 
in the equatorial region. The electric potentials measured at the surface could then directly 
reveal the super-rotation region. The numerical results of Cupal et al. PJJj suggest that 
these additional Lorentz forces could also amplify the super-rotation. 

Even when the outer sphere is at rest, we obtain slightly different results depending 
upon the direction of differential rotation. The equatorial symmetry could be broken by 
electric currents crossing the equatorial plan. This hypothesis could be tested by measuring 
the electric potential in both hemispheres in future experiments. 

Finally, we note that the ultrasonic Doppler velocity measurements show that the flow 
is not steady in the equatorial region, except when the outer sphere spins. Instead, rapid 
oscillations of the radial velocity are observed, even for the lowest |A/|. 

The latitudinal profiles of angular velocity inferred from the electric potentials are re- 
markably similar across a wide range of parameters. This seems to indicate that a rather 
universal behavior is behind all this. More numerical modelling, relaxing the assumptions 
of El and close to the experimental parameters are needed to explain our experimental 
findings. 
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We have obtained excellent results with ultrasonic Doppler velocimetry. Only radial 
velocities can be measured in the present set-up, but we are currently adapting it to ac- 
commodate measurements of the azimuthal velocities. This direct measurement will give us 
the amplitude and the radial distribution of the time-averaged angular velocity at a given 
latitude. Measurements of the induced toroidal magnetic field inside the sphere could also 
help characterize the latitudinal variation of the angular velocities. 

Global rotation appears to enhance super-rotation, in particular at large Af. As a 
consequence, very large azimuthal velocities are inferred from electric potentials, yielding 
Rm values in excess of 30, with a power input below 4 kW. This result and the observation 
of a sizeable meridional circulation make the project of building a dynamo experiment based 
on the rotating spherical Couette flow appealing |12j . 
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